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The conductance fluctuations for various types for two- 
and three-dimensional disordered systems with hard wall and 
periodic boundary conditions are studied, all the way from the 
ballistic (metallic) regime to the localized regime. It is shown 
that the universal conductance fluctuations (UCF) depend 
on the boundary conditions. The same holds for the metal to 
insulator transition. The conditions for observing the UCF 
are also given. 

PACS numbers: 71.30.+h, 71.55.Jv 

The influence of the boundary conditions (be) on 
critical phenomena in disordered mesoscopic systems 
has been demonstrated by studies of the- conductance 
distributiorJin 3 ] and energy level statistics. □ The ensem- 
ble average of the logarithm of the conductance, (\n(g)), 
is smaller for hard wall boundary conditions than for 
periodic boundary conditions. The variance (ln 2 (g)) — 
(ln(g)) 2 on the other hand is larger for hard wall than 
for periodic boundary conditions. The distribution of 
nearest neighbour energy level separations P(s) becomes 
more Wigner-Dyson-like for periodic boundary condi- 
tions. Thus, systems with periodic boundary conditions 
exhibit a "more metallic" behaviour than those with hard 
wall boundary conditions. Different boundary conditions 
also lead to different values for the universal conductance 
fluctuations (UCF) in the diffusive metallic regime where 
the mean free path I is much smaller and the localization 
length £ much bigger than the system size. The variation 
of ensemble fluctuations of the conductance as the dis- 
order increases throughout the metallic regime has been 
studiedjQ'Q but the influence of the boundary conditions 
in the ballistic regime where I exceeds the system size, 
and close to the localized regime where £ becomes com- 
parable with the system size, has not been studied in 
detail. Also, the role of the correlation length in samples 
with a true metal-insulator-transition (MIT) — i.e. in 
systems where the localization length is "infinite" in the 
metallic regime — has not been discussed. We present 
here numerical studies of cubic systems and squares with 
spin-orbit scattering — all of which have a true MIT — 
for both hard wall and periodic boundary conditions in 
the direction(s) perpendicular to transport. 

We are using the tight-binding model with the Hamil- 
tonian 
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where n, n' are nearest neighbour lattice sites on a square 
or cubic lattice. For systems with spin-orbit-interactions 
t and t' take on values of +1 or —1 and the hopping 



integrals V n ,n' ar e 2x2 matrices; without spin-orbit- 
interactions, the hopping integrals are scalar and the spin 
"variables" have only one value. We take the site energies 
e n (independent of t) to be random variables, chosen 
from an interval [—IF/2; IF/2] with a uniform probability 
distribution. The parameter IF serves thus as a measure 
of disorder strength. The conductance is calculated using 
the transfer matrix method and the Landauer formuIa.Ej 
We have used the analytical derivation of Lee et alB to 
calculate theoretical values for the UCF for both types of 
boundary conditions. In order to change the boundary 
conditions to periodic, one needs to make the following 
changes (references to equations are from the Appendix 
of Lee et al. (1987)13): in the eigenfunctions Q m to the 
diffusion equations, the cosines in the transverse direc- 
tions (Eq. (A9)) must be replaced by exponentials with 
a factor of 27r instead of tt in the argument; this will lead 
in effect to a factor 4 in the m 2 and m 2 terms in the 
modified eigenvalues A m (Eq. (A13)), and to a summa- 
tion over all integers (including negative ones) for m x and 
m y in Eqs. (A15), (A16), (A24) and (A25). The results 
are presented in Table |. The values are only half those 
given by Lee et al. due to a factor 2 in the definition of g. 
Also, our result for the three-dimensional case is slightly 
higher, probably due to our calculating the involved sums 
to a higher precision. The boundary conditions have of 
course no effect for the quasi-one-dimensional case. The 
values given here are those for the standard deviation a g 
for the orthogonal universality class of Random Matrix 
Theory. The values for the other universality classesli2l 
are obtained by dividing the variance, i.e. er 2 , by the uni- 
versality class parameter (3, where (3 — 1,2,4 for the 
orthogonal, unitary, and symplectic universality classes 
respectively. 

In Fig. m we show the standard deviation of the con- 
ductance in ensembles of 10,000 samples for different sys- 
tem sizes (squares with L x L lattice sites) and bound- 
ary conditions (open symbols: periodic boundary condi- 
tions; filled symbols: hard wall boundary conditions) as 
a function of the inverse of the average conductance. It is 
well known that in a two-dimensional disordered tightp. 
binding model all the-£jiates are exponentially localized.El 
The typical structureSfl of the fluctuations with increas- 
ing disorder strength can be seen: after an initial strong 
increase in the ballistic regime (large (g)) it reaches a 
peak value which becomes more pronounced for larger 
systems; then the fluctuations drop back to the universal 
value and finally decrease again in the strongly localized 
regime. The boundary conditions have apparently no in- 
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fluence on the behaviour outside the region of UCF in 
this case. Notice that for the case of periodic boundary 
conditions for the large system size of L = 128, a g ap- 
proaches the theoretical value of 0.393 given by the lower 
horizontal line. 

In Fig. | the same data is plotted for a two-dimensional 
disordered system of size L = 64 with periodic bound- 
ary conditions. In the same plot the mean free path and 
the localisation length as a function of disorder strength 
W are given. Both the mean free path and the localiza- 
tion length were-iobtained from the numerical results of 
Economou et alfl The localization length was obtained^ 
by the transfer matrix method, while the mean free path, 
was obtained by the coherent potential approximation^ 
(CPA). Notice that £ is always larger than I. So for a 
given system size (L — 64 in this case), there is a finite 
region where I -C L -C £■ Only in this region there is a 
plateau visible at the correct UCF-value. For W < 1,1 is 
larger than L, and we are in the ballistic regime where one 
observes a monotonic increase of a g followed by the char- 
acteristic maximum as the system enters the crossover 
between the ballistic regime and the regime character- 
ized by UCF. 

Fig. | shows that the same overall behaviour is ob- 
served also for three-dimensional systems. As there is a 
MIT (indicated by vertical lines in Fig. ||) where the con- 
ductance distribution and therefore also its standard de- 
viation become universal, i.e. independent of system-size 
(though still depending on the boundary conditionsaEl) , 
this value is approached after leaving the region of UCF. 
A direct comparison of the results shows again that the 
boundary conditions have only minimal effect outside 
that region. The additional peak noticeable in some of 
the periodic boundary conditions data are due to a near- 
degeneracy of eigen-energies for very small disorder. 

In Fig. |i| we plot the data for one of the systems again 
as a function of W, together with the mean free path I 
and the correlation length £. In the three-dimensional 
disordered case, the mean fpe path and the correlation 
length £ were again obtained!!!! by the CPA and the trans- 
fer matrix method. The size of the cube is L = 16 and 
periodic boundary conditions were used. In the three- 
dimensional case, I drops as the disorder strength W in- 
creases, while £ increases as W increases. The plateau in 
o g is seen only when I <C L and £ < L, so that we expect 
a wider plateau for larger systems. The fluctuations be- 
gin to approach the critical value as soon as L w £. In 
the three-dimensional case too, when W < 2, I is larger 
than L and the ballistic regime is observed, followed by 
the maximum in the crossover regime. 

Finally, Fig. || shows data for square systems with 
spin-orbit-intetjactions. We have chosen the Evangelou- 
Ziman model,II3 where even in the absence of diagonal 
disorder there is disorder in the hopping matrices V nin >, 
which accounts for the fact that the fluctuations do not 
vanish for small diagonal disorder. Apparently, boundary 
conditions have a noticeable influence on the fluctuations 
even outside the region of UCF, but this is likely due to 



the peculiar overall structure in this case, most signifi- 
cantly the fact that the UCF value is much smaller than 
the critical value, causing another increase in the stan- 
dard deviation as one approaches the MIT. 

In conclusion, we have investigated the conductance 
fluctuations for various types of systems with both hard 
wall and periodic boundary conditions from the ballistic 
regime to the localized regime. The boundary conditions 
seem to have a relevant influence on the conductance fluc- 
tuation only in the region of UCF and at the critical point 
of the MIT. In true metallic systems, the fluctuations 
begin to deviate from the UCF-value and approach the 
critical value as soon as the correlation length approaches 
the system size. 
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FIG. 1. Standard deviation a g of the conductance for 
squares of L x L lattice sites. Full symbols: hard wall be; 
open symbols: periodic be. Note that the be seem to have lit- 
tle effect outside the plateau region. The two horizontal lines 
indicate the theoretical values for the UCF for hard wall (top) 
and periodic (bottom) be. 
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FIG. 2. Standard deviation of a square of 64 x 64 lattice 
sites together with the mean free path I and the localisation 
length £. 
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FIG. 3. The standard deviation for systems of L x L x L 
lattice sites: a) hard wall be; b) periodic be. The horizontal 
lines indicate the UCF values; the vertical lines indicate the 
MIT. 
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FIG. 4. The standard deviation of a cube with 16 x 16 x 16 
lattice sites together with the mean free path I and the cor- 
relation length £. 




FIG. 5. The standard deviation for square systems of L x L 
lattice sites with spin-orbit-interaction (Evangelou-Ziman 
model) for hard wall (solid symbols) and periodic (open sym- 
bols) be. The two horizontal dashed lines indicate the the- 
oretical values for the UCF for hard wall (top) and periodic 
(bottom) be. 



TABLE I. The universal conductance fluctuation values for 
different be and dimensionality of the system. 



be 


[ Q1D 


2D 


3D 


hard wall 


0.365 


0.431 


0.559 


periodic 


0.365 


0.393 


0.471 
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